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ABSTRACT 

We present an analysis of the physical conditions in an extreme molecular cloud in the Antennae 
merging galaxies. This cloud has properties consistant with those required to form a globular cluster. 

We have obtained ALMA CO and 870/im observations of the Antennae galaxy system with ^ 0".5 
resolution. This cloud stands out in the data with a radius of < 24 pc and mass of > 5 x 10^ Mq. 

The cloud appears capable of forming a globular cluster, but the lack of associated thermal radio 
emission indicates that star formation has not yet altered the environment. The lack of thermal 
radio emission places the cloud in an early stage of evolution, which we expect to be short-lived 
(< 1 Myr) and thus rare. Given its mass and kinetic energy, for the cloud to be confined (as its 
appearance strongly suggests) it must be subject to an external pressure of F/ks ^ 10^ K cm“^ - 
10,000 times higher than typical interstellar pressure. This would support theories that high pressures 
are required to form globular clusters and may explain why extreme environments like the Antennae 
are preferred environments for generating such objects. Given the cloud temperature of ^ 25 K, the 
internal pressure must be dominated by non-thermal processes, most likely turbulence. We expect 
the molecular cloud to collapse and begin star formation in < 1 Myr. 

Subject headings: g3lax.ies: clusters: general; galaxies: individual(NGG 4038/9); galaxies: interactions; 
galaxies: star formation, submillimeter: galaxies 


1. INTRODUCTION 
1.1. Globular Cluster Formation 

Globular clusters are among the most ancient objects 
in the universe, often with ages >12 Gyr (Bolte & Hogan 
1995 Garretta et al.|20QQ ) and are c ommon around mas - 
sive galaxies in the universe today (Harris et al. 2013). 
The present-day abundance of globular clusters is re¬ 
markable given that that the fraction expected to survive 
10 Gyr is extremely small, potentially lower than 1% 


(iFall fc Zhang|[2QQT Whitmore, Ghandar, fc FaII||2007 ). 

Thus, this extreme type of star formation may have been 

a critical mode in the early evolution of today’s massive 
galaxies. 

Initial theories about globular cluster formation sug¬ 
gested that these objects were among t he first to gravi¬ 
tation ally collapse in the early universe (Peebles & Dicke 
1968). Subsequent work, particularly after the launch 
of the Hubble Space Telescope, has demonstrated that 
clusters with extrem e stellar densities of ten exceeding 
^ 10^ stars pc“^ (e.g. Miocchi et al.||2Q13|) can still form 
in the universe today (|0’ Connell et al. 1994) - the so- 
called “super star clusters” (SSCs). Since that time, 
numerous studies have indicated that the properties of 
SSGs are consist ent with those expected o f young globu¬ 
lar clusters (e.g. |McLaughhn & Fall||2008|). 

A number of physical processes contribute to the de¬ 
struction of clusters, including two-body relaxation, stel¬ 
lar mass loss and feedback, compressive and tidal shocks 
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as clusters orbit their host galaxy, and tidal truncation. 
The extent to which each of these processes act on a spe¬ 
cific cluster will depend on a variety of factors, including 


the o rbital properties of the cluster (Gnedin & Ostriker 
1997). Indeed, the extent to which cluster disruption is 


mass- depen dent is still debated ( Fall et al.[2QQ9[|Bastian 


et al. 2012). For unresolved clusters in galaxies outside 
the local group, there is typically limited (if any) dynami¬ 
cal information, and for clusters younger than a few Myr, 
little dynamical evolution will have taken place. For all 
of these reasons, it is not possible to say whether any 
particular SSG will survive for a Hubble time. 

While there is no generally accepted definition of “Su¬ 
per Star Gluster”, here we adopt a definition based on a 
cluster having the potential to evolve into a globular clus¬ 
ter, regardless of whether or not it actually will do so over 
the following ^10 Gyr. This requirement results in both 
mass and radius limits on the range of objects that can be 
considered as SSGs. Specifically, most present-day glob¬ 
ular clusters have half-light radii of < 10 pc (although 
some have radii as large as ^ 15 pc, van den Bergh et 


al. 1991 


and stellar masses of > 10^ Mq ([Harris et 
In addition, these clusters are expected to 
lose > 1/2 of their mass due to dy namical effects over 
10^^ years (iMcLaughlin & Fall 2008), which suggests that 
to be a globular cluster progenitor, a young star cluster 
should have a mass of > 2 x 10^ M^. If star formation 
efficiency is ^ 20 — 50% ( Ashman fc Zepf||2001 ; Kroupa 


et al. 2001), the initial molecular core from which the 
cluster is formed must have a mass of > 10^ Mq. 

Optical techniques have been able to probe the evo¬ 
lution of SSGs (and presumably some future globular 
clusters) to ages as young as a few million years. Be¬ 
fore this time, the clusters can be significantly shrouded 
by their birth material, limiting the usefulness of optical 
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observations. Beginning in the late 1990’s efforts began 
to observe SSC evolution at even earlier ages (< a few 
million years) by using radio observations to detect the 
free-free emission from the ionized gas around t he clus¬ 


ter and internal to the cluster's dust cocoon (Turner, 


et al. 

20001 |Beck, Turner, & Kovo 2000 

Johnson et al. 

200TT 

Johnson 20021 1 Johnson et al.|2003t| 

Beck et al.|2004i 


Turner fc Beck 2QQ4| [Johnson et al. 2QQ4| Reines et al 
20081 IJohnson, Hunt, fc Reines]2QQ91 Tsai et al. ]2QQ9 ' 


Aversa et al.||2011t |Kepley et al.||2014p. We refer to 
these objects as “natal” SSCs, meaning that the clus¬ 
ters themselves have already formed, but they have not 
yet emerged from their birth material. Studies of natal 
clusters were able to place constraints on the relative life¬ 
time of this enshrouded phase of SSC evolution to < a 
million years a nd the gas density of t he ionized hydrogen 
Ue > 10^ cm“^( Johnson et al. 2003). A large number of 
subsequent studies have now identified additional com¬ 
pact thermal radio sources in a number of galaxies, al¬ 
though their l ow detection rate supports their rela tively 
short lifetime ( |Tsai et ar|2009 Aversa et al.||2011 ). 

However, determining the physical conditions that give 
rise to SSCs (and their surviving descendants - globular 
clusters) has been mired in the fundamental difficulty 
that once an SSC is present in the molecular cloud, it will 
dramatically alter it. Thus in order to observationally 
probe the conditions capable of creating an SSC requires 
not only identifying mole cular clouds that are compact 
(radii < 25 pc, see Section 3.5) and massive {'> 10^ Mq), 
but also for which massive stars have not begun to dis¬ 
rupt the environment. Efforts to observe the actual for¬ 
mation of SSCs - before the star clusters have formed 
- requires high spatial resolution millimeter observations 
to determine the physical properties of the material from 
which the clusters will form. Such work has largely been 
stymied by the available observing facilities and limited 
to only the most nearby galaxies. 

One example of a relatively nearby starburst system 
in which some progress has be en made is M82. At 
3.6 Mpc (Freedman et al. 1994), relatively good linear 
resolution was achievable even before ALMA. This sys¬ 
tem was observed using the Owens Valley Radio Obser¬ 
vatory (^VRO}in^CO( 2-1) with a linear resolution of 
17 pc Keto et al. (2005). While the compact molecular 


clouds observed in 1V182 are likely to be associated with 
early SSC evolution, multi wavelength observations sug- 
gest that these clou ds have already begun star formation 
(Keto et al. 2005), and have therefore disrupted their 
birth environment^. 

This paper reports results from an ALMA Early Sci¬ 
ence project studying the Antennae galaxies. In a high 
resolution survey of C O emission from the Antennae 
(Whitmore et al. 2014) the most immediately striking 
feature was a compact, high line width cloud with lit¬ 
tle associated star formation. This is coincident with a 


sou rce ide ntified by Herrera et al. (2011); Herrera et 


al. (2012) as a potential proto-SSC using H 2 and ear¬ 


ner, lower resolution ALMA data. The strong compact 
H 2 emission appears to be due to warm (1700-2300 K) 

^ A possible exception to this is a CO cloud located at 
09h55m54.5s +6 9d40’50”, however the properties of this cloud are 
not provided by |Keto et ah] ( |2005| . 


TABLE 1 

Interferometric CO Observations of the Antennae 


Eacility 

Transition 

Beam 

Reference 

CMA 

12CO1-0 

4.91" X 3.15" 

Wilson et al. 

|2000l 

SMA 

12C03-2 

1.42" X 1.12" 

|ueda et ai.|| 

J0I2| 

SMA+PdBI 

12C02-1 

3.3" X 1.5" 

|wei et ai.|(|2 

10121 

ALMA-SV 

12C02-1 

1.68" X 0.85" 

Espada et ai. 

12012 


ALMA-SV 

12C03-2 

1.05" X 0.60" 

tlerrera et ai. 

2012 

) 

ALMA-CycO 

12C03-2 

0.56" X 0.43" 

tnis paper 


shocked gas. However, the size of the cold molecular 
component could only be constrained to < 100 pc, which 
precluded a conclusive identification. 

The present paper characterizes this source, which we 
consider among the best candidates for a proto-SSC, and 
lays out the evidence for and against the source’s even¬ 
tual evolution into a SSC or GC. 


2. ALMA OBSERVATIONS 

We obtained ALMA observations of the Antennae sys¬ 
tem in CO(3-2) and 870 /im continuum with the goal 
of probing the conditions of cluster formation and early 
evolution; data calibration is discussed in detail in the 
overview paper ( Whitmore et al.]|2014 ). Briefiy, the ob¬ 
servations consisted of a 13-point mosaic, and were car¬ 
ried out in the “extended” configuration, with a maxi¬ 
mum baseline of ^ 400 m and 5 km s“^ spectral channels. 
The resulting rms was determined using line-free chan¬ 
nels and found to be 3.3 mJy/beam. With an angular 
resolution FWHM of 0."56 x 0."43 (59 x 45 pc), these 
observations are well-matched to the expected diameter 
of the precursor giant molecul ar c louds of < 50 pc (or 
a radius of 25 pc, see Section 3.5). For this paper, we 
also recalibrated and reimaged the SV data, as well as 
compared it to previous results to check for consistency. 
The SV CO(2-l) data used in this paper was taken with 
a beam size of 1.68" x 0.85". The rms of the SV CO(2-l) 
observations was also measured using line-free channels 
and found to be 6.5 mJy/beam. 

The ALMA observations enable the study of a compact 
and luminous source in the CO(3-2) data cube (see Fig¬ 
ure . This cloud is part of the s uper giant molecula r 
cloua complex known as SGMC2 ( Wilson et al.jgOOO ); 
the full 3-D cube of SGMC2 is shown in LJgure The 
Antennae galaxy system was previously observed in both 
CO(2-1) and CO(3-2) by ALMA as part of the “science 
verification” (SV) process , which has already resulted in 
publications (see Table p!| |Espada et al.||2012t |Herrera et 

Ml2^. 

The specific source discussed here was singled out 
in the lower-resolution ALMA science verification data, 
with the CO(3-2) emission being coincident with strong 
H 2 emissio n - potentially indi cating shocks due to in¬ 
falling gas (jHerrera et al. 2012). Even with lower resolu¬ 
tion data ('^ 100 pc), it was speculated that this region 
might contain an SS C in the early stages of its evolution 
(Herrera et al. 2012). However, the spatial resolution of 
the SV CO(3-2) data is ^ 100 pc, or roughly twice that 
of the Cycle 0 data presented here. It is clear that in the 
SV data, the molecular cloud that is the subject of this 
paper is not resolved and is blended with other molecular 













































































































3 


-18®52'00.0" 


20 . 0 " 


40.0" 


53'00.0" 


20 . 0 " 






CO(3-2) 

moment 

^ I 




\ 3.6cm 
emission 


)JCO{3-2) 

emission 


Cycle 0 beam^ 


55.20s 
RA (J2000) 


12h01m51.60s 


Fig. 1.— The candidate proto super star cluster molecular cloud is located in the “Overlap” region of the Antennae merging galaxies, 
(left) HST color image of the central area of the Antennae Overlap region (red = Pact, green = F814W (~ I-band), blue = F435W (~B- 
band)), (middle) A magnified view of the “overlap” region overlaid with contours (2, 4, 8, and 16 Jy beam“^ km s“^) from the ALMA 
Cycle 0 CO(3-2) moment 0 map. (right) A zoomed-in view of the region surrounding the candidate proto super star cluster molecular 
cloud overlaid with both CO(3-2) contours of the molecular cloud after extraction from the data cube (0.4, 0.8, 1.6, 3.2 Jy beam“^ km 
s“^) and 3.6cm radio emission (-1, 4, 5, 7, lOcr, a = 3.8 X 10“^ mjy beam“^). The synthesized beam, shown in the right panel, has a size 
of 0."56 X 0."43. 



Fig. 2.— A 3-D view of the data cube illustrates the spatial and velocity structure of the entire SGMC2 region. The candidate proto-super 
star cluster molecular cloud is indicated by the arrow, and the secondary cloud is spatially offset slightly to the north-west; the primary 
cloud stands out in the data cube due to its round shape, compact size, and larger velocity dispersion (it is visibly extended in velocity 
space, in contract to other clouds in the cube, including the secondary cloud which is offset to higher velocity). The spatial dimensions of 
the cube are 9."6 x 9."6, and the velocity axis spans 1300-1800 km s“^. 


material in the vicinity. 

2.1. Cloud Analysis 

Determining the properties of this cloud requires that 
it first be isolated from other emission in the region. As 
shown in Figures and there is a redshifted secondary 
velocity component along this line of sight, that must 
be deconvolved from the primary source before analy¬ 
sis. We extract the primary cloud from the velocity cube 
for further analysis by creating a sub-cube around the 
cloud with dimensions of - 1" X 1" X 300 km s-^ The 
integrated intensity contours of the extracted primary 
cloud are shown in Figure and the observed properties 
are listed in Table We determine the cloud size by 
deconvolving the synthesized beam from the extracted 


source, which yields a half-light radius for the cloud of 
< 24 ± 3 pc (< 0.23"). The derived properties of the 
cloud are given in Table As a sanity check, the prop¬ 
erties of the primary cloud were also determined using 
the CPROPS program on the entire data cube (no t ex¬ 
clusively the sub-cube) ( [Rosolowsky fc Leroy||2006 ), and 
the resulting parameters agree to within the uncertain¬ 
ties. Thus, by-hand measurement of the half-light size, 
automated Gaussian fits, and moment based measure¬ 
ments all yield roughly consistent sizes for our cloud. As 
this is a marginally resolved object with a clearly mea¬ 
surable line width methodological uncertainties do not 
overwhelm any of our conclusions. Throughout this text, 
we refer to the properties of the extracted primary cloud 
only, unless noted otherwise. 
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TABLE 2 

Measured Properties of the Molecular Cloud 


R.A. 

(J2000) 

Dec. 

(J2000) 

(km s-i) 

Sc03-2 
(Jy km s“^) 

av 

(km s-l) 

Size 

FWHM (arcsec^) 

12:01:54.73 

-18:52:53.2 

1524T3 

52T5 

49T3 

0.66 ±0.12 X 0.55 ±0.07 


Note. — Properties o f the clouds were measured both with the CPROPS program 
|Rosolowsky & Leroy|(|2006|), and with Gaussian fitting (FWHM = 2.35cr). The quoted uncer- 
tamties retlect empirically determined variations in these values for different fitting attempts. 
The synthesized beam has not been deconvolved in this table to strictly report measured 
properties. The deconvolved size is reported in Table 



Fig. 3.— The CO(3-2) and CO(l-O) spectra of a - 1" x 1" re¬ 
gion around the candidate proto super star cluster molecular cloud 
taken from the full data cube (i.e. the primary cloud has not been 
extracted). The line emission indicates that there are two com¬ 
ponents along the line-of-sight that have different temperatures. 
Spectral profile of the CO(3-2) line from current work (ALMA 
Cycle 0 obs ervations) and the C O(2-l) line from ALMA science 
verification jEspada et al.||2012|). The CO(3-2) data is convolved 
to the CO(2-i) beam and both data sets are corrected for beam 
dilution. There is clearly more than one velocity component; in 
these convolved data, the secondary source was fit by a Gaussian 
and subtracted from the spectra. The CO(3-2)/CO(2-l) ratio is 
dramatically different in the two velocity components. 


2.2. Relative Astrometrie Solutions 

A comparison between the CO(3-2) emission and data 
at other wavelengths, requires an understanding of the 
relative astrometric accuracy. Based on the phase stabil¬ 
ity of the ALMA observations, we estimate the absolute 
astrometric accuracy to be better than ^ 0.2". Cen¬ 
timeter observations from the VLA have an astrom etric 
accuracy better than ^ 0.1" (Brogan et al. 2010), and 


therefore the 3.6 cm and CO(3-2) observations have a 
relative precision of better than the synthesized beam 
of the ALMA data, and we consider them to be astro- 
metrically matched. We also register the astrometry of 
archival Hubble Space Telescope observations shown in 
Figure by matching the Paa emission throughout the 
Antennae system to common features in the 3.6 cm emis¬ 
sion. 


3. RESULTS 

3.1. Cloud Temperature and Optieal Depth 

We constrain the temperature and optical depth of the 
cloud using the CO(3-2) and CO(2-l) emission. We re¬ 


trieved, recalibrated, and reimaged CO(2-1) observations 
from ALMA’s science verification period, shown over¬ 
plotted in FigureThe CO(3-2) observations were con¬ 
volved to the synthesized beam of the CO(2-1) observa¬ 
tions and corrected for beam dilution using the Cycle 0 
CO(3-2) source size, resulting in peak brightness temper¬ 
atures of T 3_2 = 17 ± 3 K and T 2-1 = 18 ± 3 K. The 
largest angular scale to which the CO(3-2) observations 
are sensitive is ^ 6", and therefore we do not expect 
that any flux is resolved out on the size scales of interest 
here. The secondary component has a significantly lower 
CO(3-2)/CO(2-l) ratio, indicating much cooler gas than 
the primary cloud. __ 

RADEX non-LTE modeling (Ivan der Tak et al.|2007) 
was used to analyze the CO(3-2) and CO(2-l) emission. 
The line intensities and their ratio were compared to a 
grid of RADEX models covering a range of values for ki¬ 
netic temperature, CO column density, and H 2 volume 
density. The best-fit values result from a chi-squared 
minimization. Since the source is marginally resolved in 
CO(3-2), for CO(3-2) we set the beam filling fraction to 
1, and for CO(2-1) we set it to the dilution factor, or the 
ratio of the CO(3-2) size to the CO(2-l) beamsize. The 
RADEX models indicate that these transitions are opti¬ 
cally thick - the best fitting depth is r ^ 3.5T0.5, but 
the data do not rule out significantly higher values. The 
lines appear to be close to thermalized, with an excita¬ 
tion temperature within a degree of the kinetic temper¬ 
ature of 25^2^ K; this temperature is on the upper end 
of the range of th ose found for dense molecular clouds 
in the Milky Way (Shirley et al. 2013). However, there 
is a degeneracy between the inferred temperature and 
density of the cloud, and the cloud could be warmer for 
densities ^ 6 x 10^ cm“^. In other words, the observed 
brightness could also be reproduced with a large column 
of subthermally excited, warm, relatively diffuse gas. 

The temperature inferred here for the CO cloud is 
roughly 100 x less than that inferred for the compa ct H 2 


emission observe d in this region of 1,700-2,300 K ( Her 
et al. |20TT ). In addition, the H 2 EW HM line-width 

1^0 U O ~ 1 r-r/^T*T*/^T*0 / 001 1 \ TO oirv-mi-R 


(2011) is significantly 


rera _ _ 

of ^ 150 km s“^ Herrera et al. 
higher than the FWHM line-width of the CO emission 
measured here of 115 km s“^ (c^ = 49 km s“^). There¬ 
fore we infer that the origin of the CO and H 2 emission 
may not be identical. We speculate that the H 2 emission 
has a low filling factor, sampling only the most strongly 
shocked regions in and/or around the cloud. 


3.2. Cloud Mass 
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TABLE 3 

Derived Properties of the Molecular Cloud 


Deconv. 
FWHM (pc2) 

Half Light 

R (pc2) 

^^virial 

(10® M 0 ) 

Mxco 

(10® Mo) 

M.jiADEX 
(10® Mo) 

^Cont. 

(10® Mo) 

Kin 

(K) 

P 

(g cm-3) 

< 53 X 41 

< 27 X 21^ 

29-85 

3.3 - 15 

> 2.8 

3.1 - 7.4 

23-35 

2 - 20 X 10-21 


Note. — The FWHM as measured after deconvolving the synthesized beam and adopting a distance of 
21.5 Mpc. The half light radius resulting from deconvolution assuming a Gaussian distribution of light. The 
cloud is consistent with either being marginally resolved or a point source at this resolution. 
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Fig. 4.— These observations allow us to disentangle the 
molecular cloud associated with the secondary velocity compo¬ 
nent. Contours of the CO(3-2) moment 0 map of the molecu¬ 
lar cloud extracted from the 3D data cube (0.4, 0.8, 1.6, 3.2 Jy 
beam“^ km s“^) overlaid on the color moment 0 map of the data 
cube with the primary cloud extracted. The cloud corresponding to 
the secondary velocity component can be seen in spatial projection 
with the primary cloud. 


The mass of the cloud is estimated using four different 
methods, each subject to different caveats. First, given 
the source size, velocity dispersion of = 49±6 km s“^, 
and assuming an isothermal sphere we calculate the virial 
mass to be in the range of = 2.9 — 8.5 x 10^ Mq, 
which is consistent with the virial mass estimated from 


ALM A SV observations of ^ 5 x 10^ Mq (Herrera et al 


2012). This mass estimate will only be valid if the cloud 


IS m gravitational virial equilibrium; any additional ve 
locity in the cloud will result in this method overestimat¬ 
ing the mass. Given the complex dynamical structure of 
SGMCs, we treat the estimated virial mass as an upper 
limit. 

For the second method we use the RADEX models of 
the CO(3-2) and GO(2-l) observations to determine the 
best fit column density of N^o ^ 5 x 10^^ cm~^. If we 


assume an abundance ratio of nH^I'^c o = 2 x 10^ (Blake 
et al. II 1987 Wilson & Matteucci |1992 ), this column den- 


sity of CO corresponds to an H 2 mass of y^non-LTE = 
2.8 X 10^ Mq. However, the x^ values are shallow toward 
higher values of Nqo, and only weakly constrain the up¬ 
per bound. The models also indicate that the CO(3-2) 
has an optical depth of r > 3, and therefore this mass 
estimate is a lower limit. 


The third method we employ assumes a conversion 
factor, Xqo- Values for Xqo la starbu rsts are known 


to va ry by at least a factor of four (Bolatto et al 


2013), and thus the mass estimated using Xqo should 


be regarded as uncertain by a corresponding factor 
Here we adopt a “starburst” CO-to-H 2 con version fac¬ 
tor X^g = 0.5 X 10^^ cm“^(K km s“^)“^ (Bolatto et 


al. 2013). If we were to adopt a “standard” Xqo = 


2 X 10^ cm“^(K km the resulting cloud mass 

would be a factor of four larger. We further assume that 
CO(3-2) is thermalized with respect to CO(1-0) given the 
brightness temperatures of T 3_2 = 17±3 K and T 2-1 = 
18 ± 3 K ; if the lines are not thermalized and CO(3-2) 
is relatively underpopulated, this mass estimate will be 
low. This method results in Mxco — — 15 x 10^ Mq. 

The last method estimates the dust mass from the 



the range of mass values determined above, we adopt a 
cloud mass of M= 0.3 — 1.5 x 10^ Mq. We note that the 
virial mass that would be inferred for this cloud appears 
to be a factor of 5-10 x too high. 

Given the estimated mass and size of this cloud, the 
resulting gas volume density is p > 100 Mq pc“^. While 
there are molecular clouds found in the Milky Way with 
masses of ^ 10^ Mq, their radii are 2 — 4x larger, re¬ 
sulting in significantly lower surface densities. Similarly 
massive clouds have also been identified in other nearby 
galaxies ( Bolatto et al.|2008 Donovan Meyer et al.|: 


their surface densities are also far lower (see h'igure 
We estimate the mass of the stellar cluster that will 
potentially result from this molecular core by assum¬ 
ing a star formation efficiency (SFE, fraction of mass 
turned into stars over the lifetime of a cloud). The 
net efficiency can vary wildly, ranging from a few per¬ 
cent in the Milky Way (e.g. |Lada & Lada| |2003|), to 
> 50% in cluster-forming cores by basic boundedness 
arguments. If we adopt a SEE typical for clusters of 

s) of ^ 20% —50% (jAshman & 


Zepf 

2001 

Kroupa et al. 2001), an initial cloud mass of 

5 X ] 


.0^ Mq would result in a cluster with a mass 


of Mstar = 1 — 5 X 10^ Mq. This would be among the 
most massive SSCs that have formed in the Ante nnae if 
it forms a single cluster (Whitmore et al. 2010). Even 
if the star formation efficiency were as low as ^ 5%, the 
resulting cluster would have a mass > 2 x 10^ Mq, still 
in the regime of super star cluster masses. 
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Fig. 5. — Detecting the continuum emission associated with this 
source allows us to determine its dust mass. Contours of the CO(3- 
2) moment 0 map of the primary cloud overlaid on a color-scale 
image of the 870 p,m continuum. 

3.3. Constraining the Ionizing Flux Associated with the 

Cloud 

In order to assess the extent to which star formation 
may have already affected the physical state of the ion¬ 
ized gas, we searched for ionizing flux potentially coming 
from stars within the molecular cloud. Figure shows 
the Paa emission in the region, and while there is diffuse 
emission associated with the SGMC in general, there is 
no discrete source associated with the molecular cloud 
in question. However, it is possible that given the em¬ 
bedded nature of this source, Paa emission could suffer 
from significant extinction, and thus we also utilize radio 
observations to identify free-free emission that might be 
present. 

We created two radio maps from archival 3.6 cm VLA 
observations (proposal codes AN079, AP478, AS796, 
AA301); one with a synthesized beam of 0.65" x 0.42" 
to best match the beam of the ALMA CO(3-2) observa¬ 
tions and resolve out diffuse emission, and a second with 
a synthesized beam of 1.12" x 0.85", which has greater 
sensitivity to extended emission. In the higher resolu¬ 
tion map, there is no discrete source coincident with 
the molecular cloud discussed here, and the 5cr detec¬ 
tion threshold of the radio emission corresponds to an 
ionizing flux of N Lyc ^ 6 x 10^^ s~^. Th is is equivalent 
to ^ 60 0-type stars (Vacca et al. 1996). This limit is 


roughly three times lower than the previous limit placed 


on possible ionizing flux in this region o f 2 x 10^^ s ^ 


using lower resolution 6 cm observations (Herrera et al. 
2011|). For comparison, the ionizing flux o f 30 Dor has 

4 X 10^^ s“^ (Crowther & 


3een estimate d to be ^Lyc 
Dessart|[l998 ). 


In the lower resolution map, while there is no discrete 
source coincident with the molecular cloud, there is dif¬ 
fuse emission in the region. Without velocity information 
for this diffuse emission, it is not possible to disentangle 
potential line-of-sight confusion in this complex region. 
However, to estimate the possible contribution from stars 
within this molecular cloud to the diffuse emission, we 
first fit and subtract Gaussian profiles to the other dom¬ 
inate sources in SGMC2. After subtracting the emission 
due to nearby sources, we measure the flux density due 
to diffuse emission in the cloud aperture to be 0.03 mJy, 


which corresponds to an ionizing flux of 1. 1 x 10^^ s~^, 
which corresponds to ^ 100 0-type stars (Vacca et al. 
1996|). However, the diffuse morphology of this emission 


IS not consistent with it coming from a single compact 
source. 

Thus, while both the Paa and cm radio observations 
show diffuse emission associated with the SGMC2 struc¬ 
ture, there is no discrete source associated with the com¬ 
pact cloud discussed here (Figure!^. We constrain the 
possible ionizing flux that could oe due to embedded 
stars in this cloud to be < 60 0-type stars. We can 


> 10^ MO ( 

Leitherer et al. 

1999 

), roughly two orders 

of magnitude smaller than t 

le am 

ticipated cluster mass 


(Section 3.2). Either star formation has not begun or it 
is so deeply embedded that its ionizing radiation is con¬ 
fined by gas continuing to accrete onto the protostars. 

3.4. Determination of Cloud Pressure 

With a surface density of ^ 4 x 10^ Mq pc“^ and size- 
linewidth coefficient of cr^/i? = 90 km^ s“^ pc, the cloud 
is not consistent with being in either pure gra vitational 
virial eq uilibrium or free-fall collapse (Figure [^|Heyer et| 
al.||2009|). Nevertheless, the morphology of the cloud in- 
^ates that gravity is playing a significant role (round, 
compact, and bright - making it stand out as a singu¬ 
lar object in the data cube), suggesting that it is not a 
transient object. As illustrated in Figurethe observed 
line width value can be explained if the cloud is subject 
to external pressures of F/ks ^ 10^ K cm“^, roughly 
five orders of magnitude higher th an that typical in th e 
interstellar medium of a galaxy (Jenkins et al. 1983). 
This is consistent with theoretical considerations that 
have argued SSC form ation requires extrerne pressures 


i^/ks ^ 10^ K cm ^) dJog fc Soloman[|l992| Elmegreen 


fc EfremLOv|1997 Ashm^fc Zept|2QQip . This high pres¬ 
sure is also in accord with previo us hndings of co mpres- 
sive shocks in the overlap region ( Wei et ^|2Q12 ) 


3.5. Expected Proto Super Star Cluster Cloud Size 

The expected physical size of a molecular cloud capable 
of forming an SSC can be estimated based o n virial the¬ 
orem arguments. Eollowing previous work (jElmegreen 
1989), the external cloud pressure Pe, cloud mass A/, 
cloud radius r, and velocity dispersion Gy can be related 
by 

mual 

47rr3 ’ 

where H is defined by = n(ne), and here we adopt H = 
0.5. It has been estimated that SSC formati on requires 
internal pressure s of Po/k > 10^ K cm“^ (Elmegreen 
& Efremov 1997). These high pressures inhibit the dis¬ 
persal of the natal material and achieve sufficiently high 
star formation efficiencies in the cloud core. If we adopt 
a minimum mass of M = 10^ Mq for a cloud capable of 
forming a SSC, the resulting cloud radius is r ^ 25 pc. 
Likewise, for the velocity dispersion of Gy = 49 km s“^ 
observed for the cloud discussed here and the apparent 
external pressure of P^/k > 10^ K cm“^ (see Eigure 
the expected radius is r ^ 25 pc, which is within the 
uncertainty of the cloud half-light radius of < 24 ± 3 pc 
extracted from these observations. 

4. DISCUSSION 
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Fig. 6. — The measured properties of the primary molecular cloud iudicate that it may be subject to pressures of /ks = 10^ K cm“^. 
This plot shows the relatiouship betweeu surface deusity aud the size-liue width coefhcieut with the caudidate proto-super star cluster cloud 
plotted as the blue poiut (error bars represeut Icr). Liues correspoudiug to pure gravitatioual virial equilibrium (uo exterual pressure) aloug 
with correspoudiug hue for equilibrium with exterual pressures of P/A:^ = 10®, 10^ , 10^, aud 10^ K cm~^ are showu iu black, correspoudiug 
couditious for clouds uudergoiug free-fall collapse are showu as red dashed liues [Field, Blackmau, &: Ketol (|2011|). For comparisou, data 
poiuts are showu from mole cular cloud survey iu uearby galaxies as light greeu poiuts (|Douovau Meyer et al.||2ui3|, aud Galactic clouds 
are showu as purple poiuts ( |Heyer et al.|200^ . 


4.1. On the Origin of the High Pressure Inferred for the 
Moleeular Cloud 

If the cloud is confined, as we expect from such a strong 
concentration of gas, then external pressure must play a 
key role. This high external pressure could result from 
the weight of the surrounding SGMC2 (if the SGMC is 
roughly in hydrostatic equilibrium) and/or large scale 
compressive shocks. The pressure generated by the 
surrounding molecular material can be estimat ed using 
Pg/kn ~ 1-5 c in~3 K (McioMd/Mo)2(r/pc)“^ ( |Bertoldi 
& McKee 1992). The SGMC2 region has an estim ated 
total mass of ^ 4 x 10^ Mq (Wilson et al. 2000) and 
radius of ^ 400 pc, resulting in a pressure from the over- 
lying molecular material of F/ks ^ 10^ K cm“^, which 
is at least an order of magnitude less than the internal 
pressure inferred for this cloud. 

Given that this cloud is not only in the “overlap” re¬ 
gion of the Antennae, but also appears to be at the nexus 
of two filaments of GO(3-2) emission that are suggestive 
of colliding flows (Whitmore et al. 2014), a significant 
amount of external pressure could also be generated by 
ram pressure. This interpretation is consistent with pre¬ 
vious work indicating strong H 2 line emissi on and an 
abrupt ve l ocity gradient acros s this region (Herrera et 


al. 


2011 


Herrera et al. 


2012). However, the morphol- 


ogy of the cloud suggests an isotropic source of pressure, 
which is in tension with a ram pressure origin. 


4.2. Timeseales for Cloud Evolution and Star 
Formation 

The relevant timescales for this cloud to evolve are 
driven by the free-fall time and the crossing time. The 
compact size and marginally-resolved round morphology 
suggest that self-gravity has had a significant role in 
shaping the source, although that is not possible to con¬ 
clusively demonstrate with the data in hand. The cloud 
must be largely supported by turbulence: given the in¬ 
ferred density of the cloud (n^ 10^ cm“^), if the pressure 
were entirely thermal, it would require a gas temperature 
of ^ 10^ K, which is not reasonable or consistent with 
these observations that indicate T ^ 25 K. We conclude 
that this cloud is most likely supported by turbulence. 
This conclusion is also supported by the virial mass be¬ 
ing nearly an order of magnitude larger than the mass 
inferred from the dust continuum, indicating significant 
internal motion contributing to the observed line width. 

Thus, we adopt the turbulent crossing time as the ap¬ 
propriate timescale for the evolution of this cloud and 
estimate it as ter ^ Djay ^ 1 Myr, where D is the di¬ 
ameter of the region. If self-gravity is not important, the 
cloud will disperse on this timescale - the turbulence will 
dissipate on this timescale in any case. If self-gravity is 
important (as argued above), the cloud will collapse on 
this timescale. The free-fall time of the cloud can be 
estimated as t/j = {Sir/32 G ~ 8 x 10^ years. 

Given the strong associated H 2 emission (Herrera et al. 






























Johnson et al. 


2012) and line-width, which indicates internal velocities 
higher than virial equilibrium, it is plausible that this 
cloud has already begun free-fall collapse, in which case 
we are witnessing a very short-lived stage of cluster evo¬ 
lution. Given these arguments, we expect this cloud to 
collapse on timescales < 1 Myr. 

4.3. Expected Number of Proto-SSC Molecular Clouds 

Given that this phase of SSC formation is expected to 
be extremely short-lived (< 10^ — 10^ yr), it is not sur¬ 
prising that this is the only clear example that we have 
found to date of a compact cloud without associated star 
formation that is sufficiently massive to form an SSC. 
Within ^ 1 Myr, we expect that this cloud will be asso¬ 
ciated with star formation, similar to the clouds observed 


in M82 by Keto et al. (2005). 

The estimated current star formatio n rate fSFR) for 
the Antennae is ^7 — 20 yr~ ^ (Zhang, Fall, & 
Whitmore 2001 Brandi et al.j2009). if cluster forma¬ 
tion follows a power -law distribution of dN/dM oc M~^ 


(Zhang & Fallll999) with lower and upper masses of 10^ 


and 10^ Mq, we expect ^ 20% of the stellar mass to 
be formed in clusters with M> 10^ Mq. Thus we ex¬ 
pect only a few SSCs masses of > 10^ Mq to form every 
^ 5x 10^ years. Optical studies indicate that there are six 
SSCs with ages < 10^ years and masses > 10^ Mq (Whit- 
more et al.|2010 ). This suggests that a massive cluster is 
formed every <10^/6 ~ 1.7x 10^ years, in agreement with 
the estimate derived above from the total SFR. The pre¬ 
dicted number of pre-stellar molecular clouds capable of 
forming an SSC with mass > 10^ Mq can be generalized 


- /SFR(Meyr-i)\ 

^VSSC-GMC — 0.2 X tcollapse I -- I Cj 

where_Coiiapse is the timescale for cloud collapse (see Sec¬ 
tion 4.2). Given a dynamical timescale for this cloud of 


^ 0.5 — 1 Myr, we expect that if we observed the An¬ 
tennae system at any point in its recent star forming 
history 10^“^ yr), we would find at most one cloud of 
this evolutionary state and mass. 


4.4. Implications for Globular Cluster Formation 


The physical properties of this cloud can provide in¬ 
sight into a mode of star formation that may have been 
dominant in the earlier universe, when globular clusters 
were formed prolifically. In particular, the pressure of 
P/ks ^ 10^ K cm“^ supports the hypothes is that such 


high pressures are necessar y to form SSCs (Ashman & 
Zepf||2QQT Elmegre^|2QQ2 ). Based on what appears to 


be a nearly universal power-law distribution of cluster 
masses, numerous studies have suggested that the for¬ 
mation of SSCs is statistical in nature, resulting from a 
“size of sample” - galaxies with higher star formation 
rates will form more clusters overall, and the formation 
of SSCs r esults from populating the tail of the mass dis¬ 
tribution ( Fall fc Chandar|2Q12 ). Based on these ALMA 
observations, we suggest an alternate interpretation: in 
a turbulent interstellar medium, the pressure distribu¬ 
tion also has a power-law form, and the properties of 
clusters that form track the pressure distribution. Thus, 
while galaxies present a power-law distribution of clus¬ 
ter masses, only regions with sufficiently high pressure 
will be able to form the most massive SSCs. If pressures 
P/ks ^ 10^ K cm“^ are indeed required to form SSCs 
(and the surviving globular clusters), similarly high pres¬ 
sures must have been common during the peak of glob¬ 
ular cluster formation a few billion years after the Big 
Bang. 
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